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1 Introduction 



Functional equations methods appeared in the theory of exactly solvable lattice models 
intimately connected with Baxter's commuting transfer matrix method [lj. In the early 
seventies Baxter introduced in his pioneer work pQ the concept of Q-operators and T — Q 
relations determining the eigenvalues of the transfer matrix of the corresponding vertex 
model. 

The complex calculations involved in Baxter's construction of Q-operators seems to 
have restricted its use and other functional methods, such as the Reshetikhin's analytical 
Bethe ansatz, were employed instead to obtain the spectrum of transfer matrices related 
with quantum Kac-Moody algebras [2j. However we remark here the recent progresses in 
the construction of Q-operators by employing quantum algebras representation theoretic 
methods [HI S]. 

On the other hand, the advent of the algebraic Bethe ansatz in the late seventies pro- 
vided a systematic approach to find the eigenvalues and eigenvectors of transfer matrices of 
integrable vertex models. This method, originally proposed by Takhtadzhan and Faddeev 
[5], is based on the existence of a pseudovacuum or reference state and appropriate commu- 
tation rules provenient from the Yang-Baxter algebra, which is a common algebraic structure 
associated with integrable vertex models. 

Although the algebraic Bethe ansatz method is a powerful tool exhibiting a rich math- 
ematical structure, its implementation for models that do not possess a trivial reference state 
is still an obstacle to be overcomed. The aim of this paper is to show that the Yang-Baxter 
algebra can be explored in order to generate functional relations determining the spectrum of 
transfer matrices. In order to illustrate that we consider the XXZ model with non-diagonal 
twists and open boundaries for general values of the anisotropy and boundary parameters. 
Besides the relevance of studying models with general boundary conditions in the context 
of statistical mechanics [6], the XX Z model with non-diagonal twists and open boundaries 
are both included in a class of models where a trivial reference state is absent. 

This paper is organized as follows. In section 2 we describe the XXZ model with 
general toroidal boundary conditions. In particular, we discuss the case of non-diagonal 
twists and we derive functional relations for its eigenvalues making use of the Yang-Baxter 
algebra. In section 3 we approach the eigenvalue problem for the XXZ model with non- 



1 



diagonal open boundaries using the algebraic-functional method devised in the section 2. 
Concluding remarks are discussed in the section 4 and in the appendices A-D we give some 
extra results and technical details. 

2 The XX Z model with non-diagonal twisted bound- 
ary conditions 

The advent of the Quantum Inverse Scattering Method [51 Ej was an important stage in the 
development of the theory of exactly solvable quantum systems. This method unveiled a deep 
connection between solutions of the Yang-Baxter equation, quantum integrable systems and 
exactly solvable lattice models of statistical mechanics in two dimensions [S]. In statistical 
mechanics an important role is played by vertex models whose respective transfer matrix is 
constructed from local Boltzmann weights contained in a operator C A j. Let V be a finite 
dimensional linear space, the integrability of the vertex model is achieved when the operator 
valued function £ : C — > End(V (g> V) is a solution of the Yang-Baxter equation, namely 

£i 2 (A - n)C 1 z{\)C 2 M = £ 2 3<»£i3(A)£i2(A - //), (1) 

defined in the space V\ ® V 2 <g> V 3 . Here we use the standard notation dj e End(Vi ® Vj). 
The complex valued operator C A j(X) can be viewed as a matrix in the space of states A 
denoting for instance the horizontal degrees of freedom of a square lattice, while its matrix 

L 

elements are operators acting non-trivially in the j-th position of VJ. In its turn the space 

i=l 

Vj represents the space of states of the vertical degrees of freedom at the j-ih site of a chain 
of lengh L. 

The transfer matrix of the corresponding vertex model can be written in terms of the 
monodromy matrix Ta(A) defined by the following ordered product 

T A {\) = £ A l{\)Cal-i(\) ■ ..C A1 (X). (2) 

As a consequence of the Yang-Baxter equation, the monodromy matrix satisfies the following 
quadratic relation 

R(X - fj.) T A (X) ® T a (jjl) = T A (/j,) ® T A (X) R(X - fi), (3) 
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usually denominated Yang-Baxter algebra. The i?-matrix appearing in ([3]) follows from the 
solution of the Yang-Baxter equation through the relation -R(A) = P£(A) where P denotes 
the usual permutation operator. The above f?-matrix plays the role of structure constant 
for the Yang-Baxter algebra and it consist of an invertible complex valued matrix acting on 
the tensor product A® A. 

The invariance of the Yang-Baxter algebra plays an important role in the description 
of integrable spin chains with general toroidal boundary conditions. One can easily verify 
the invariance of ([3]) under the tranformation T A {X) — > G A T A {\) provided that the c-number 
matrix G A is a symmetry of the .R-matrix, i.e. 

[R{\),G A ®G A ] = 0. (4) 

Consequently we can define the operator 

T(X)=Tt a [G a T a (X)] (5) 

which constitutes an one parameter family of commuting transfer matrices, i.e. [T(A), T(/i)] = 
0. 

Now restricting ourselves to the XXZ model, V{ = C 2 and the corresponding L- 
operator is that of the anisotropic six vertex model 



£{\) 



fa(\) 








\ 





6(A) 


c(A) 








c(A) 


6(A) 





V o 








a(\)J 



(6) 



whose Boltzmann weights are given by a(A) = sinh(A+7), 6(A) = sinh(A) and c(A) = sinh(7). 

In the Ref. [UJ UH] the authors discuss the possible classes of twist matrices G A com- 
patible with the i?-matrix associated with (J6j). These twist matrices turn out to be 



G A =\ {it) G A = (7) 





where a and (3 are arbitrary complex parameters. The above matrices G A are non-singular 
for a, (3 7^ and the XXZ model with generalized toroidal boundary conditions is obtained 
as the logarithmic derivative of the transfer matrix §5§ at the point A = [HI ITU] . 

The understanding of physical properties of vertex models demands the exact diago- 
nalization of their respective transfer matrices, which can provide us information about the 



free energy behaviour and the nature of the elementary excitations. When the twist matrix 
of type (i) is considered, a trivial reference state is available and the eigenvalues of the cor- 
responding transfer matrix can be obtained through the algebraic Bethe ansatz j5j [7j. By 
way of contrast, the twist matrix of type (ii) breaks the U(l) symmetry of the the system 
leaving only a Z 2 invar iance. 

Although a trivial reference state in the framework of the algebraic Bethe ansatz is 
no longer available when the matrix G_a of type (ii) is considered, the eigenvalues of the 
corresponding transfer matrix were obtained in [9] by means of the Baxter's T — Q method. 
It is worthwhile to remark here that the isotropic case associated with the XXX spin chain 
admits any 2x2 twist matrix and interesting enough the algebraic Bethe ansatz solution 
can be obtained by exploring the GL(2) symmetry [TT] . 

From the historical point of view, the T—Q method was introduced in Baxter's remark- 
able works on the eight vertex model [1] and more recently it has found applications in many 
areas such as the study of integrable systems [12], conformal field theory [31 IB], correlations 
functions [H] and efficient description of finite temperature properties [15J. Motivated by 
the ideas of Baxter's T — Q method [H IB] and the algebraic Bethe ansatz [SI E], in what 
follows we shall demonstrate how we can use the Yang-Baxter algebra to obtain functional 
relations determining the spectrum of the transfer matrix T(X). 

Let us consider the eigenvalue problem for the transfer matrix 



taking into account the type (ii) G_a matrix given in (JTj). As shown in the appendix A, we 
can set a — (3 — 1 without loss of generality for the purposes of this paper. Now considering 
the definition ([2]), the monodromy matrix Xa(X) consist of a 2 x 2 matrix whose elements 
are operators that we denote 



T(X) M = A(A) M 



(8) 




(9) 



Therefore, the transfer matrix §5§ reads 



T(X) 



B(X)+C(X). 



(10) 





(11) 
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is not an eigenstate of the transfer matrix (TTOT) . Nevertheless, the state |0) is still of great 
utility. Considering the definition together with (jSJ), the elements of 'Ta(X) satisfy the 
following relations 

A(A)|0) = a(A) L |0> D(X) |0) = b(X) L |0) 

B(X)\0) = f C(A)|0) = 0, (12) 

where the symbol f stands for a non-null value. Notice the relations (TT2l) imply in 

T(X) |0> = B(X) |0> (13) 

and 

T(X)B(X) |0> = B(X)B(X) |0) + C(X)B(X) |0> . (14) 
Now we have reached a point of fundamental importance. The term C(X)B(X) |0) 



present in the right hand side of (fl4l) can be evaluated with the help of the Yang-Baxter 
algebra ([3]). In order to show that we collect the following commutation rule among the ones 
encoded in the relation Q, 

C(A)fl(Ai) = B(ji)C(X) + [A(jm)D(X) - A(X)D(ji)) . (15) 

The commutation rule ffTol) . together with the relations f|T2|) . results in the following identity 

C(X)B(ji) |0) = C ^—^ [a^) L b{X) L - a(X)H(ri L ] |0) . (16) 

Therefore, as we are interested in C(X)B(X) |0), the next step is to consider the limit /i — > X 
in the above relation which can be evaluated using L'Hopital's rule. Thus we are left with 

C(X)B(X) |0> = M(A) |0) (17) 

where M(A) = Lc(A) 2 a(A) L - 1 6(A) L - 1 . 

At this stage we have already gathered the basic ingredients to obtain functional rela- 
tions determining the eigenvalues of T(A). Operating with the dual eigenvector (tfj\ on the 
left side of Eqs. ( ITBl and ( 1141) we are left with the relations 

A(A)F = F 1 (X) (18) 
A(A)Fi(A) = F 2 (A) +M(X)F , (19) 
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where F = <</>|0>, F 1 (X) = (^| B(X) |0> and F 2 (A) = {if>\ B(X)B(X) |0). 

Here we remark one of the roles played by integrability in this approach. Since the 
transfer matrix T(A) belongs to a commuting family, the eigenvectors are independent 
of the spectral parameter A. Thus the dependence of the functions i*i(A) with A is solely 
determined by the operator B(X). In the appendix C we have computed that dependence 
making use of Eqs. © and (JHJ). 

From the Eq. (1181) we can see that the functions A(A) and -Fi(A) differ only by the 
constant factor F , thus they possess the same zeroes. Now we can eliminate A(A) from 
the Eq. (|18|) and substitute the result into the Eq. (TIT?]) which yields the following relation 
involving only the functions i^(A), 

F 1 (X) 2 = F 2 (X)F + M(X)Fl (20) 

Considering the Eqs. (j2J), (JSJ) and (jHJ), in the appendix C we have demonstrated that 
the functions -Fi(A) and F 2 {X) can be written as 



m) = m) ^<E^A and m) = Fm Y^^, (2!) 

1=1 



sinh^) " Li sinh(Af )x 



where Xj denote the zeroes of the function F,(A). Now a closer look in the Eqs. (TIB"]) and 
( fT9|) reveals that we can obtain the eigenvalue A(A) by determining the variables A^ together 
with the ratios ^r^. 

Fo 

In order to obtain analogues of Bethe ansatz equations determining the variables Xj, 
we first consider the Eq. ( 120|) at the points A = A^- . Then we find 

F 2 (A; i) ) + M(a5 1) )F = 0. (22) 

Next we consider the points A = X^ in Eq. (I20p which yields 

F 1 (Af) 2 -M(Af )F 2 = 0. (23) 

Furthermore, we observe by setting A = in the Eqs. (j!8p and (I19p that the ratios can 
be written in terms of A(0) which can be easily evaluated. Therefore we are left with the 
following expression for the ratios ^W^, 



f 
F 2 (0) 



F 



A(0) 

A(0) 2 . (24) 



o 
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In order to avoid an overcrowded section we have performed the diagonalization of T(0) in 
the appendix B. The eigenvalue A(0) turns out to be 



A(0) = sinh(7) i e^ 2: 



r = 0, ...,2L- 1. 



(25) 



Gathering our results so far, from the Eqs. ( 1181122112311241) we obtain the following 
expression for the transfer matrix eigenvalues A(A), 



A(A) = smh(j) L e^ if " A) 



i=i 



sinh(A 



0-h 



(26) 



recalling that r = 0, . . . , 2L — 1 and provided that the variables A * satisfy the following 
system of non-linear algebraic equations 



smh(A^- + 7) sinh(A^ 
sinh(7) sinh(7) 



(l)^ 



L-1 



sinh(A? } - \P) 

—r- n 



i=l 



sinh(A 



(2h 



smh{\f + 7) sinh(Af j ) 
sinh(7) sinh(7) 



(2> 



L-1 



2iirr T 1 



L 



n 



sinh(Af } - Xf) 



sinh(A 



j = l,...,L-l 

(27) 

j = l,...,2(L-l) 
(28) 



Now we shall examine some aspects of the Eqs. (1271) and (12 8p taking into account the 
crossing properties discussed in the appendix D. Considering the Eq. ( 1D.9I) . it follows that 



A'(A) = (-1) X+1 A(-A - 7) 



(29) 



where A* (A) denotes the eigenvalue of the transposed transfer matrix T*(A). Although the 
relation ( 1291) does not imply in A(A) = (— l) i+1 A(— A — 7), we have verified the existence 
of eigenvalues satisfying that relation through the direct diagonalization of T(A) for small 
chain lengh L. 

Let us suppose the relation A(A) = (— l) i+1 A(— A — 7) holds for some eigenvalue. 
Hence, from the Eq. ( TPS!) we can conclude that the function Fi(A) also satisfy 



F 1 (A) = (-l) z+1 F a (-A-7) 



(30) 



Since the function M(A) enjoys the property M(A) = M{— A — 7), from the Eq. (fTQl 
we also find that 

F 2 (A) =F 2 (-A- 7 ). (31) 



Assuming that the relations (130]) and fl3T|) hold they have remarkable implications 
concerning the solutions of Eqs. (j2*7|) and Let us analyze first the odd L case. Thus 

L — 1 is an even number and the Eq. ( l30i) implies, for instance, in the following relation for 
the variables A^- , 

A^ = -A«- 7 < = l,...,^i. (32) 

Thus altogether we have only independent variables A^. 

Now we turn our attention to the L even case. In that case the Eq. ( l30i) implies for 
instance in the relation 

A«^ = -A«- 7 i = (33) 

while the remaining root X^-i * s nxe d at the value — 2. 

(2) 

The above analysis can also be carried out for the roots A] taking into account the 
Eq. ( 13TT) . Since 2(L — 1) is an even number, the variables A^ must be related for instance 
by 

AgLx = -Al 2) - 7 < = 1,...,L-1. (34) 

These implications of the crossing symmetry have been verified by a numerical analysis 
of the Eqs. ( |271) and ( 1281) for small values of L. By comparing the eigenvalues given by the Eq. 
( |2U|) with the direct diagonalization of T(X), we have also verified that the Eqs. (I26|I27|I28I) 
indeed describes a complete spectrum. 

We close this section by remarking the dependence with -j- in the right hand side of 
Eqs. (I2T1) and (!28l) . As far as we know, this kind of dependence with the chain lengh had 
not appeared previously in the context of Bethe ansatz, which may be of relevance for the 
description of the thermodynamical limit L — > 00. 

3 The XX Z model with non-diagonal open boundaries 

In Sklyanin's pioneer work [T7j, the author has been able to generalize the Quantum Inverse 
Scattering Method to accomodate integrable spin chains with open boundaries. The Yang- 
Baxter equation is still the corner stone of Sklyanin's approach and it turns out that the 
XXZ model with general open boundary conditions can be obtained from the following 
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double-row transfer matrix 



t(u) = 



Tr A \K\(u)T A (u)K A (u)T A (u) 



(35) 



where T A {u) = C AL (u)C AL -i{u) . . . C A \(u) and f A {u) = C Al {u)C A2 (u) . . . C AL (u) are the 
standard monodromy matrices that generate the corresponding closed spin chain with L 
sites. The integrability at the boundaries is governed by the matrices K A (u) and K A (u), 
each one describing the reflection at one of the ends of an open chain. 

Moreover, the boundary conditions compatible with the bulk integrability are con- 
strained by the so-called reflection equations, which for K A (u) reads 



C 2 i(u - v)K 2 (u)£ 12 (u + v)K x (v) = K x (v)C 2 i{u + v)K 2 {u)C 12 (u - v), 



(36) 



while a similar equation should also hold for the matrix K A (u). 

Turning our attention to the XXZ model, the bulk hamiltonian is described by the 
^-operator (JH]) and the most general boundary matrices satisfying the reflection equations 
are of the form 

k 2i(u) k 22 (u) j A \k 21 (u) k^ 2 (u) 

whose matrix elements are given by 



K A {u) 



(37) 



kn( u ) = smh(hi + u) 

k 2i( u ) = ^3 sinh(2«) 

k xl (u) = sinh(/i^ — u — 7) 

k ti( u ) = sinh(— 2u — 27) 



ki 2 (u) = h 2 sinh(2ii) 
k 22 {u) = sinh(/ij" — u) 
k X2 {u) = h 2 sinh(— 2u — 27) 



k 22 {u) = sinh(/i^ + u + 7). 



(3? 



The above K- matrices possess altogether six free boundary parameters {hf} and we find 
the XXZ model hamiltonian with open boundaries, 

L-l 



K = E *X+i + a i a i+i + c osh(7Kaf +1 



4 = 1 



+ 



sinh(7) 
sinh(ft,J~) 
sinh(7) 
sinh(/i^") 

related to the former by 



(h 2 + h 3 )al + i{h 2 - h 3 )a{ + cosh(/i 1 )a{ 
04 + ht)al + i(ht - h+)al + oosh(fc+K 



(39) 



t\0) = 2 sinh(/^) smh(hf ) sinh(7) 2L_1 cosh(7)ft + L cosh(7) 2 + sinh(7) 



(40) 



where af, of and erf denote the usual Pauli matrices acting on the z-th site. 

For general values of the boundary parameters hf the double-row transfer matrix 
does not exhibit £7(1) invariance, which makes the application of the algebraic Bethe ansatz 
quite complex due to the lack of a trivial reference state. By imposing hf = hf = the 
matrices K A (X) (|37p become diagonal and the £7(1) symmetry exhibited by the bulk of the 
system is restored. In that case the spectrum of (1351) was obtained in [T7] through the alge- 
braic Bethe ansatz. However, several progresses concerning non-diagonal open boundaries 
have been reported in the literature over the past years when the boundary parameters or 
the bulk anisotropy satisfy a given constraint [T%]-|31]. 

In what follows we shall explore the Yang-Baxter algebra to generate functional rela- 
tions determining the complete spectrum of the double-row transfer matrix ( l35i) for general 
values of the bulk anisotropy 7 and boundary parameters hf. In order to tackle the eigen- 
value problem 

t(u) |$) = A(u) |$) (41) 

for the double-row transfer matrix (1351) . we first notice the ^-operator ([6]) satisfies the prop- 
erty C\2{u)C\2{— u) = a(w)a(— w)Id<8> Id, where Id denotes the 2x2 identity matrix. There- 
fore, the monodromy matrices T A (u) and l~ A {u) are related by T A [u) = [a(u)a(— u)] L 7^ 1 (— u) 
and we find the following relation, 

f A {u) ® Id] R(2u) [T A {u) <g> Id] = [Id ® T A (u)) R(2u) [id ® T A (u)] , (42) 

with a straightforward manipulation of the Yang-Baxter algebra ([3]). Analogously to Eq. 
©, we denote the matrix elements of T A (u) by 

( A(u) B(u) 

t A (u) = I . (43) 

1 ' \C{u) D(u) ' 

In this way the double-row transfer matrix ( J35l) reads 



t(u) = kf^U^u) + k+(u)U 2 (u) + kf 2 (u)U 4 (u) + kf^ujU^u) (44) 

where 

Uxiu) = k n (u)A(u)A(u) + k u (u)A(u)C(u) + k 21 (u)B(u)A(u) + k 22 (u)B(u)C(u) 
U 2 (u) = ki X (u)C(u)B(u) + ku(u)C(u)D(u) + k 21 {u)D{u)B(u) + k 22 (u)D(u)D(u) 
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U 3 (u) = k n (u)A(u)B(u) + k u (u)A(u)D(u) + k 21 (u)B(u)B(u) + k 22 (u)B(u)D(u) 
U 4 (u) = k n (u)C(u)A(u) + k u (u)C(u)C(u) + k 21 (u)D(u)A(u) + k 22 (u)D(u)C{u). 

(45) 

Next we consider the action of t(u) on the state |0) defined in (TTT]i . which will require the 
following relations 



A(u)\0) = a(u) L \0) D(u) |0> = b(u) L |0) 
B{u) |0) = f C(u) |0) = 



(46) 



obtained from the structure of the monodromy matrix T^(u), and the relations (1T2"]) from 
the previous section. 

Then, in order to evaluate t(u) |0), we shall make use of the commutation rules between 
the elements of T^iu) and 7a(u) enclosed in the relation fj4"2l . In particular, we shall use 
the following ones, 



A(u)B(u 

D(u)B(u 

A(u)B(u 

C(u)B(u 
B(u)B(u 



K2u) -B(u)A(u) - C -^B(u)D(u) 



a(2u) 



a(2u) 



a(2u) -, . , c(2u) . „. , 
^B{u)D{ U ) + ^A{u)B{u) 

b{2u) B(u)A(u)- C -^B(u)D(u) 



a(2u) 
c(2u) 



a(2u) 
B(u)B(u). 



a(2u) 

A(u)A(u) - D(u)D(u)} + B(u)C(u) 



(47) 

Considering the Eqs. fj4*4"|) . fT45l) . ffT2"]) and fT4oT) together with the relations fj4"T|) . a straight- 
forward calculation leave us with 



t{u) |0> = Y {u) |0> + Y l {u)B{u) |0> + y x («)B(u) |0> + Y 2 {u)B{u)B{u) |0> 
where the functions Y (u), Yi(u), Yi(u) and Y 2 (u) are given by 



(46 



a 2L (« 



— b 2L {u) sinh(u + 7 + /i^) sinh(-u + 7 — /i x 



sinh(/i 1 l_ — u) sinh(/i 1 + w) — 



sinh(2u + 27) 



sinh(2-u + 7) 

sinh(2n) 



sinh(2w + 7) 

a L (u)b L (u) (h 2 h^ + h 3 hf\ sinh(2w) sinh(2w + 27) 



Y 1 (u) 



a ( u 



)h 3 sinh(/i+ — u) sinh(2 



^sinh^w + 27) 
sinh(2« + 7) 
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_r, \,i . i / , s . , ln s sinh(2w + 27) 

sinh(2u + 27) 



Y\{u) = — a L (u)h^ sinh(/i 1 + «) sinh(2-u) 



sinh(2« + 7) 
sinh(2w + 27) 



+ b (u)h 3 sinh(w + 7 + sinh(2u) 
K 2 (m) = -/4/iJsinh(2w)sinh(2w + 2 7 ). (49) 

We can now follow the scheme devised in the previous section and operating with the dual 
eigenvector ($| on the left side of Eq. PH|) we obtain 



A 



[u)f = Y Q (u)f + Y 1 (u)f 1 (u) + Y^Uu) + Y 2 (u)f 2 (u), (50) 

where f Q = ($|0>, h{u) = ($| B(u) |0>, h{u) = ($| |0> and / 2 («) = ($| B(u)B(u) |0). 

Before we proceed it is important to examine first the Eq. f[5*U]) under the light of the 
results obtained in the appendix D. Indeed, the crossing symmetry enables us to estabilish 
the relation B(u) = (—l) L+l B(—u — 7) ( ID. 81) which implies in 

f 1 (u) = {-l) L+l f 1 {-u- 1 ). (51) 

Furthermore, the commutation rule B(u), B(u) = given in (HTl) together with the relation 
(ID. 81) . imply in the crossing invariance of the function f 2 {u), i.e. 

f2(u) = f 2 (-u-j). (52) 

From the results of the appendix D we also have the crossing invariance of the double-row 
transfer matrix t(u) = t(—u — 7), which implies in 

A(u) = A(-u-7). (53) 

In their turn the functions Yi(u) given in f)49p fulfill the relations Y (u) = Y (—u — 7), 
Yi(u) = (— l) L+1 li(— u — 7) and Y 2 (u) = Y 2 (—u — 7). Taking into account the above 
properties, one can easily verify the crossing invariance of the Eq. ( 150]) . which simplifies to 
the following functional relation, 

A(u)/ = Y {u)f + Y l (u)f 1 (u) + Y 1 (-u - 7 )A(-m - 7) + Y 2 {u)f 2 {u). (54) 



Considering now the results from the appendix C, the function fi(u) (1C 13[) can be 
written as 

j- 1 \ i' /a\ ^vr sinh(uf ) -u) 

fx{u) = fi(0) I] . ■ (55) 

i=i smh(i^ ') 
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On the other hand, the functions f 2 {u) and A (it) given by (IC.14I) and (IC.16I) will receive 
extra simplifications due to the relations (|52|) and (|53|) . Thus they can be written as 

(56) 
(57) 



, . , • i / \ • , / x tt sinh(u| - u) smh(u\ + 7 + u) 

f 2 (u) = g 2 smh(w) smh(« + 7) [[ — - (2) 7TT1Z) — 

smh(-Uj- ) smh(ii- +7) 



8=1 



/\( \ a fn\ TT sinn (^i 0) - u ) sinh(?4 0) + 7 + n) 

A(«) = A(0) 1| ; . w (0) , x 



i=i 



The initial condition A(0) has been determined in the appendix B through a direct 
analysis of t(0). Then the next step is to determine the variables . In order to do that 
we consider the Eq. (|5"4l at the points it = it = Uj and it = uf\ By doing so we find 
the following system of algebraic equations, 



n(«f)/o + ri(<0/iW^ = (58) 

YoiuV) - A(uf ] )} f + Y 1 (-uf ) - ^M-uf - 7) + Y 2 {uf)h{uf) = (59) 

fi{-uf-l) = 0.(60) 

A direct inspection of the Eqs. fl58H60p shows that the ratios ^j^- and ^ are also required. 
In contrast to the case of twisted boundary conditions considered in the previous section, 
where M(0) = 0, the point u = of Eq. (|54|) is not very enlightening. However, in order to 
determine the ratios and ^ we consider, for instance, the non-trivial points u — U\ and 
u = u 2 such that Y\{u\) = Y\(u 2 ) = 0. Then we are left with the following equations 



,(o> 



,(0)x 



,(°) 



(0) 



,(0)x 



Y Q {uf) - A(uf) / + Y x {uf)h{uf) + Yi{-uf - i)h{-uf - 7) 



,(2) 



(2)^ 



.(2) 



/i(0) 



\Y (ui) - A(ui)] f + Y 1 {-u i - 7 )/i(-«i - 7) + y 2 («i)/ 2 (Mi) = 



1,2 



(61) 



which can be solved for the required ratios. 

Now our results can be gathered and we are left with the following expression for the 
eigenvalues A (it), 



2L ^ji 2 sinh(iij- Uj — it) Bvah.{u^ } + 7 + u 



,(°) 



A (it) = 2 cosh(7) smh(hf) sinh(/t 1 ) sinh(7) 2i JJ 



(o) 



i=l 







sinh(u 



(Oh 



sinh(U'°' ) + 7) 



, (62) 



provided that the variables u%, u 2 \^ and uf 1 satisfy the following system of non-linear alge- 



braic equations, 



sinh(itj + ^) 



sinh (ui — 



h 3 sinh {ui — ^ — hf 
hs sinh(i/j + | — K{ 



i = 1,2 



(63) 



1 We have performed the shifts m —> m — \. 
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y&P) n 



L - x sinh^ - uf ] ] 



>i(-^ 0) -7)n 



L-l 



i=l 



sinh(u 

^(u^') sinh(wy U ' ) ) sinh(w 



• 1/ (i) . . (o^ 



(o) 



7) II —712) 

i=i sinn(u) 



i=i 

.„<•»; 



sinh(w. 



det(ifi) 



sinh(4 2) + 7 + w$ 0) )det(#2) 



sinh(wf ) + 7) det(Ho) 

j = l,...,L + 2 



(64) 



2 cosh( 7 ) sinh(tf) sinhM sinh( 7 ) 2i if " ^ + 7 + ^ - 



sinh(wf') sinh(-Uj-°' ) + 7) 



i=l 

-Y 2 (u) ') sinh(u} smh( W j + 7) U . ~ ( (2) ' . w (aT7 r Hptr/y n 

sinn(u> ; ) sinh(«; ; + 7) Q.et{n ) 



i=l 

a, A ^ sinh( ^ (1) + j + uf ] ) det(ffi) 
i=i sinh(^ 1) ) det(^o) 



(65) 



o v,^ • • ufh-\ ■ U ^2L ft sinh (^ 0) - U T) sinh(^ U; + 7 + uf) (2 -, 
2cosh(7j smn(/ii) sinh(/i 1 ) smh(7) 7^ ^ Yq{u ' 



i=l 



sinh(u^) sinh(-Uj-°' ) + 7) 



^ sinh^-^f) (2) / rT 1 sinhK (1) + 7 + ^ 2) ) 

i=i smh(^ ; ) i= i sinn(it- 'J 



det(ifi) 
det(tfo) 

j = l,...,L-2 

(66) 

where are 2x2 matrices resulting from the solution of the system of equations ( 16T1) for 
the ratios and a r. These matrices turn out to be 

JO JO 

( (j)l{ux) 2 (wi)\ ( 4>o{ui) <j>2{ui)\ ( 4>l{ui) 0o(wi)\ 

#0 = i?! = # 2 = (67) 

\0i(m 2 ) h(u 2 ) J V ^0(^2) <fo(uz) ) \<M« 2 ) 0o(wa)/ 

where 

, / x , / s • -, /, +n • 1 /, s • 1 / tt sinh(^ 0) + 2 - it) sinh(^ 0) + 2 + it) 
Vu) = 2 cosh( 7 ) smh(/i+) sinh^ ) smh( 7 ) 2L n \ \ f 2 ^ (0) 2 



i=i sinh(ii' ) sinh(w- + 7) 
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j L ^ sinh^^ + i + 



/ i=i sinh(w> J 

I / \ w 7\ • i / 7 U . u 7^ 2 s inh K (2) + f -M)sinh( M f ) + ^ + M 
&2(w) = l2(«- -)smh(w- -)smh(M+ -) [[ 



2 2 2 j=1 sinh(M J - 2 ' ) ) sinh(M 4 - 2 ' ) + 7) 

(68) 

Taking into account the Eq. ( HU1) we are left with the following expression for the 
eigenergies E of the hamiltonian 



L+2 1 sinhf'-yl 2 

E = -sinh( 7 ) 2 £ . (ok ■ (o)- , " ^osh( 7 ) - (69) 
i=i sinn(«- ; ) smn(«- ; + 7) cosn^j 

given in terms of the roots u\°\ We end this section by remarking that numerical checks 
performed for L = 2,3,4 show that the spectrum generated by the Eqs. fl62l) - fl66l) is complete. 



4 Concluding Remarks 

In this paper we have proposed a functional method in the theory of exactly solvable models 
based on the Yang-Baxter algebra. Using this method we were able to derive the eigenvalues 
of the XX Z model with non-diagonal twists and open boundaries for general values of the 
bulk anisotropy and boundary parameters. Our solution is presented in terms of analogues 
of Bethe ansatz equations whose variables involved are precisely the roots of the transfer 
matrix eigenvalues and roots of auxiliary functions defined in terms of the monodromy 
matrix elements. 

Concerning the XX Z model with non-diagonal twists discussed in the section 2, we 
stress the unusual dependence of the Eqs. ( )2*7j) and (12*8|) with the chain lengh L. We hope 
the computation of physical properties in the thermodynamical limit L — > 00, such as the 
interfacial tension obtained in [9J, will be benefited by the use of Eqs. (12 7[ I2B1 . As shown in 
|34j . we also remark that our solution corresponds to the eigenvalues of Baxter's eight vertex 
model with elliptic modulus k — 1 and a particular choice of diagonal twist matrix. 

Regarding the case of non-diagonal open boundaries described in the section 3, we 
remark the solution recently presented for general values of the anisotropy and boundary 
parameters obtained from the representation theory of g-Onsager algebra [33] and the multi- 
ple reference state structure found in [321 EH] • An interesting problem would be unveiling the 



15 



connection between these approaches. Since our method is based on the Yang-Baxter alge- 
bra, which is a common algebraic structure underlying integrable vertex models, we expect 
that our approach can be applied for other models with general open boundary conditions 
based on g-deformed Lie algebras [36J. 

Finally, we observe that our solutions involve more than one kind of variable resembling 
the so-called nested Bethe ansatz equations, which are typical of models with higher rank 
symmetry We remark that a similar result had been reported previously in the literature 
obtained from generalized T — Q relations [27] . 
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Appendix A: The condition a = (3 = 1 

Let T(X) denote the transfer matrix 

f(A) = Tr A [g a C al (X)C al ^(X) . ..C A1 (X)} (A.l) 
where £^(A) is the /^-operator given in (El) and G A is the following twist matrix 

^[m]- (a ' 2) 

(yfr \ 

Now considering the matrix M. = one can verify that M. A G A M. A = ya/3 



v^, 
where G A is given by 

o)' <A ' 3) 

which corresponds to G A with a = f3 = 1. By noticing that A / lJ 1 Ai A 1 C A j(X)Ai A Aij = 
C A j(X), we then find the relation 

M J X f W IlM^yfip T(X), (A.4) 
j=i i=j 

where T(A) is the transfer matrix fflUj) considered in the section 2. Thus the eigenvalues 
A (A) of the transfer matrix T(A) are related to the ones of T(A) by 

A(A) = v/o^ A(A). (A.5) 



Appendix B: The Eigenvalues A(0) and A(0) 

In this appendix we consider the diagonalization of the transfer matrices T(0) and t(0) 
associated with the XX Z model with twisted and open boundary conditions respectively. 

• Twisted Boundary Conditions: 

The ^-operator (EJ) consist of a regular solution of the Yang-Baxter equation, i.e. £(0) = 
sinh(7)P, where P denotes the permutation operator. In this way the transfer matrix (1101) 
at the point A = can be written as T(0) = sinh(7) L where the operator O is given by 

6 = Tr^ [GaPalPal-i ■ ■ ■ Pai] ■ (B.l) 
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Using the permutator algebra, P 2 = Id® Id and PAjPAi — PaiPiJi we can readily obtain the 
following expression for O, 

6 = G l P 1L ...P Vi P l2 , (B.2) 

and its n times product 

n 

6 H = G X G 2 . . . G n I] PjLPjL-l ■ ■ ■ Pjn+1- (B.3) 
3=1 



From the Eq. flB.3|) we have 

L = G 1 G 2 ...G L , (B.4) 
/0 1\ . . 

which implies, with G = , in 2L = I where / is the identity matrix. Thus the 

\1 Oj 

eigenvalues O of the operator O are given by 

= e^ r = 0,1,..., 2L- 1 (B.5) 

leaving us with the following expression for A(0), 

A(0) = sinh( 7 ) L e^ r = 0, 1, . . . , 2L - 1 . (B.6) 

• Open Boundary Conditions: 

The matrix K^(u) presented in Eqs. (I57|) and (1551) consist of a regular solution of the 
reflection equation, i.e. K^(0) = sinh(/ij~)Id. Together with the regularity property of the 
/^-operator ([6]) and the permutator property P 2 = Id <g> Id, the double-row transfer matrix 
( I35p at the point u = is given by 

t(0) = sinh( 7 ) 2i sinh(/i^)Tr^ [iTj(O)] /. (B.7) 

Therefore, we can see that the transfer matrix t(0) is proportional to the identity matrix I 
with the following eigenvalue 

A(0) = 2 cosh( 7 ) sinh(/i+) sinh(/i^) sinh(7) 2L . (B.8) 



Appendix C: The functions i^(A) and fi(X) 

In this appendix we derive the form of the functions Fi{\) and fi(\) used in the sections 2 
and 3, considering the algebraic properties of the monodromy matrices Xa(X) and Ta(A). 
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In order to determine the dependence of the functions i^(A) with the spectral parameter 
A, let us first introduce the following notation for the monodromy matrix ([2]) 

(L) (A L {\) B L (X)\ 

rj L) (A)= , (c.i) 

\Cl(X) D l {X)J 

differing from the one used in ([HD by the index L inserted to emphasize we are considering 
the ordered product of L matrices C_Aj{\). For instance, the matrices £.4,- (A) are given by 



where 



7i(A) 



'a(X) \ / 0' 

aj(X) = ( (3j(X) = 

" 6(A) I U(A) 0, 



j 



c(A)\ /6(A) 

^ = 1.0 o), WH (o .(A) 1 (C ' 3) 



are matrices acting non-trivially in the j-th site of the chain. 

Moreover, the monodromy matrix T^ L \X) defined by Eq. ([2D satisfies the recurrence 
relation 

Ti L+1) (A)=^ L+1 (A)Tj L) (A), (C.4) 
corresponding to the following relations for its matrix elements 

A L+1 (X) = a L+1 (X)A L (X) + (3 L+1 (X)C L (X) 

B L+1 (X) = a L+1 (X)B L (X) + (3 L+1 (X)D L (X) 

C L+ i(A) = 7l+i(X)A l (X)+5 l+1 (X)C l (X) 

D l+ i(X) = 7L+1 (X)B L (X) + 5 L+1 (X)D L (X). (C.5) 

At this stage, it is important to keep in mind that the spectral parameter A enters in 
the monodromy matrix only through the Boltzmann weights a(A) and 6(A), since c(A) is A 
independent ([6]). For the forthcoming discussion, it is convenient to introduce the weight 
u>(A) only to denote the Boltzmann weights possessing dependence with A. Thus, w(X) can 
represent both a(A) and 6(A). The reason for the introduction of w(X) will become clear 
in what follows. From the Eqs. (J2J), ( 1C.2I) and ( 1C.3I) it is clear that -Bi(A) and Ci(A) are 
polynomials of degree in w, while ^i(A) and -Di(A) are polynomials of degree 1 in w. These 
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initial conditions allow us to use the Eqs. (IC.5I) to determine the degree of the elements of 
T^\\) in w for arbitrary lengh L. Thus we obtain that -Bl(A) and C^(A) are polynomials 
of degree L — 1 in w, while the operators Al(X) and -Dl(A) are of degree L. 

In terms of the variables x = e x and q = e 7 the Boltzmann weights a(A) and b(X) 
can be written as a(A) = | (xq — x~ l q~ Y ) and 6(A) = ~ (x — x -1 ). Thus u>(A) is of the form 
-w(A) ~ x + a; -1 , and expanding its powers we find 

B{\) ~ x L - x + . . . + 

(C6) 



cl-i(x 2 ) l 1 + c l (x 2 ) l + . . . + co 



Now considering the functions Fi(A) = (V>| 5(A) |0) and F 2 (A) = (^| B{\)B{\) |0) used in 
section 2, and the fact that is A independent, we can conclude that 

Fi(A) = x-^Pxix) 

F 2 (A) = ^-^2^), (C.7) 

where Pi(x) and i^^) are polynomials in the variable x 2 of degrees L — 1 and 2(L — 1) 
respectively. Therefore, the polynomials Pi(x) and P2{x) can be written as 

L-l 



p 1 (x) = Pl n [x 2 - (xfY 

i=l 

2(L-1) 

JM*) = p 2 n [^ 2 -(^ (2) ) 2 ] (c. 



i=l 



where pj are constants and xp denote the zeroes of the polynomial Pj{x). Using the variables 



_ we th en have 



wo = w» n ( ,;,, 2 ; ' (c.9) 

i=l 



sinh(A^ 



where the constant pj is absorved by F/(0). 

Let us now consider the monodromy matrix (A) = £_ai(A)£4 2 (A) . . ■C^l(X), with 
matrix elements 

_ (L) (A L {\) B L {\)\ 

rj L) (A)= _ - , (cio) 

\Cl(X) D l {X)J 

appearing in the construction of the double- row transfer matrix ( |35i) . Compared to the nota- 
tion previously used ( )43i) . we have also included the index L stressing we are considering the 
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ordered product of L matrices Cj L j{X). The monodromy matrix T^\X) obeys the recurrence 
formula 

Ti L+1) (A)=Tj L) (A)^ L+1 (A), (C.11) 
corresponding to the following recursion relations for the matrix elements 

A L+1 (X) = A L (X)a L+1 (X) + S L (A) 7 L+i(A) 

B L+1 (X) = A L (X)P L+1 (X)+B L (X)5 L+1 (X) 

C L+1 (X) = C L {X)a L+1 (X) + D L {X)j L+1 {X) 

D L+1 (X) = C L (X)p L+1 (X) + D L (X)6 L+1 (X). (C.12) 

The same arguments previously used for the elements of (A) can now be consid- 
ered for the elements of T^\X). Thus we obtain that the operators B L (X) and Cz,(A) are 
polynomials of degree L — l in w while A L (X) and D L (X) are of degree L. Hence the function 
fi(u) = ($| B(u) |0) required in the section 3 can be written as 

h(u)=mn s ™ H f- u) . (c.i3) 

i= i smh{u\ ') 

On the other hand the function fz(u) = ($| B(u)B(u) |0) requires an extra analysis. Noticing 
that aj(0) = ^y7i(0)/?j(0) and ^(0) = (0)7,(0) together with /3|(A) = 7?(A) = 0, one 
can easily see from the Eqs. ( 1C.5I) and (1C.12I) that B(0)B(0) = 0. This same analysis can 
be performed for the point A = —7 yielding B(— ~f)B(— 7) = 0. Thus the function f2{u) can 
be written as 

2(L-2) .,/(2) \ 

h(u) = <7 2 smh(w)smh(« + 7) ]J \ , 2 , ; . (C.14) 

i=i smh(^ 

In the section 3 we also require the function A (it) defined as 

AM-figJP- (C ' 15) 

The above results together with the Eqs. (pHI) . (Ho]) , and (IHHj) show that the double-row 
transfer matrix consist of a polynomial of degree 2(L + 2) in u>. Thus the eigenvalues A(w) 
are of the form 



sinh(i4 — u) 

i=l 



a(«)=a(o) n \ ; (oh , (c-16) 

smh(ii- ) 



where m,- "* denote the zeroes of A(u) 
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Appendix D: Crossing symmetry 

In this appendix we deduce some important properties used in the sections 2 and 3 which 
are consequence of the crossing symmetry. 

Besides the Yang-Baxter equation, the £-operator §§§ satisfies other important rela- 
tions, namely 

Temporal invariance: /^^(A) = L\i (A) 

Crossing symmetry: £12 (A) = —ViC^—X — ^Vf 1 (D.l) 

( ° l \ 

where V = and U stands for the transposition in the space with index i. 

In order to estabilish a relation between the matrices T^A) and 734(A), we first recall 
the definitions 

T A (X) = £ AL (X)£ AL _ 1 (X)...£ A1 (X) (D.2) 
f A {\) = £ A1 (X)...£ AL ^(X)£ AL (X), (D.3) 

and observe that £% (A) = C% (A) which follows from the temporal invariance f ID.lh . Next 
we consider the transposition in the auxiliary space A, 

Tj? (A) = C^(X)...C%(X) 

= C%{\)...C%{\), (D.4) 

which allows us to use the crossing symmetry ( ID. II) to obtain 

TX A {X) = (-i)Lyx 1 £ M (-\-'y)V A ...VZ 1 £ AL (-\-<y)V A 

= (-l) L V^T A (-X-j)V A . (D.5) 

Now we also consider the transposition t — t\ . . Al. Then we find 

rX A \x) = [c AL (X)...c A1 (x)} tAh - tL 

= £^(X)...£^(X) 

= C Al (X)...C AL {X), (D.6) 
which leads immediately to the relation 

TX A \X) = f A {X). (D.7) 
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In terms of the matrix elements, the relations (ID.5I) and (ID .70 read 



A(X) 
B(X) 
A\X) 
B\X) 



■l) £ 5(-A- 7 ) 



l) L+1 B(-X- 7 ) 



A(X) 
(7(A) 



C(A) = (-l) i+1 C(-A- 7 ) 
D(X) = (-1) L A(-X- 1 ) 
C 4 (A) = 5(A) 
D\\) = D(\). 



(D.8) 



A straighforward calculation taking into account the relations flD.8j) reveals that the 
twisted transfer matrix ffTUl) satisfies the relation 



T*(A) = (-l) L+1 T(-A- 7 ) 



while we have the following identity for the double-row transfer matrix ([35 



(D.9) 



t{u) = t(—u — 7). 



(D.10) 
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